Universal behaviour of the 7*7 — > (77°, rj, 77/) transition form factors 
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The photon transition form factors of n, n and rj' are discussed in view of recent measurements. It 
is shown that the exact axial anomaly sum rule allows a precise comparison of all three form factors 
at high-Q 2 independent of the different structures and distribution amplitudes of the participating 
pseudoscalar mesons. We conclude: (i) The 7T7 form factor reported by Belle is in excellent agreement 
with the non-strange 1 — component of the r/ and rf form factors obtained from the BaBar 
measurements, (ii) Within errors, the n~f form factor from Belle is compatible with the asymptotic 
pQCD behavior, similar to the n and rf form factors from BaBar. Still, the best fits to the data 
sets of 7T7, 777, and 7/7 form factors favor a universal small logarithmic rise Q 2 F Pl (Q 2 ) ~ log(Q 2 ). 
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1. INTRODUCTION AND RESULTS 

The extensive experimental study of the 7*7 — > P re- 
actions (P = 7T°, 77, 77') (Tj-Q attracted much attention 
from theorists (for recent references see The 
reason was the report of the BaBar collaboration [4] 
about a persistent rise of the combination Q 2 F. K1 (Q 2 ) 
in the Q 2 region from 10 GeV 2 to 40 GeV 2 . The mea- 
sured form factor sur pas sed the asymptotic behaviour 
Q 2 F 7IJ (Q 2 ) -> v^/tt IH, U = 0.131 GeV, predicted 
by perturbative QCD (pQCD). Several theoretical inves- 
tigations 0, [IH HH, [H| indicated, however, that the cor- 
responding increase of the irj form factor for large Q 2 
values is hard to explain. Very recently new experimen- 
tal information came from the Belle collaboration [l8j]. 
The data for Q 2 F 7ri (Q 2 ) presented by this group show 
only a very mild (if any) increase in the high-Q 2 region. 

In this situation a comparison of the 7r7 form factor 
with the ones for r\ and rj can be helpful. At first sight, 
this appears difficult because of the different structure 
(different quark distribution amplitudes) of these parti- 
cles. In the present note we show however, that by using 
the exact anomaly sum rule, the high-Q 2 behavior of tt, 
r\ and 77' form factors is essentially determined by the 
well-known lowest order spectral representation of the 
triangle quark diagram. Therefore, a comparison of the 
three form factors is possible and can be trusted. 

The application of the anomaly sum rule requires to 
relate the full-QCD spectral densities with the spectral 
densities obtained from perturbative QCD. This can be 
achieved by using the concept of duality. This way the 
anomaly sum rule offers the interesting possibility 
to calculate the transition form factors without referring 
to the QCD factorization theorem. No assumptions are 
needed about the light-cone distribution amplitudes of 
pseudoscalar mesons with their specific end-point behav- 
ior and Gegenbauer coefficients. 

The high-Q 2 behavior is determined by the high-energy 
dependence of the spectral density in the corresponding 
integrals, see Eqs. (fl"U)) and ([TT]) below. At high energy — 



above the resonance region — the spectral density can be 
very well approximated by perturbation theory. This 
spectral density is the same for the three form factors. 
(The effect of the difference between the current masses 
of strange and non strange quarks should be negligible 
at high energy). Thus we can conclude that at high-Q 2 
the functional dependence of these form factors should 
be the same. A detailed analysis in [3] suggests that the 
universality may be expected already at Q 2 > 10 GeV 2 . 
Different decay constants provide for different multipli- 
cation factors but do not affect the slope. 

Looking now at the data, the following observations 
can be made: 

(i) The BaBar measurements of the 777 and 7/7 form 
factors (Fig. [1} are within errors compatible with pQCD 
factorization which implies saturation of the combination 
Q 2 Fp 7 (Q 2 ). Still, the data seem to indicate a very mild 
(e.g. logarithmic) rise with Q 2 . 

(ii) The large-Q 2 behaviour of the form factor F 7T1 (Q 2 ) 
as observed by BaBar (Fig. [5]) is in some conflict with the 
saturation predicted by QCD factorization. These data 
suggest an increase of Q 2 F 7TJ (Q 2 ). The rise as seen from 
the high-Q 2 points of the BaBar data is much larger than 
the rise observed for the 77 and 77' form factors. 

(iii) The Belle [l8[ measurement of the irj form fac- 
tor Fig. [5] is within errors compatible with pQCD fac- 
torization which implies saturation of the combination 
Q 2 Fp 7 (Q 2 ). Still, also these data seem to indicate a log- 
arithmic rise with Q 2 , but a very mild one like in the 77 
and 77' data. 

We now apply our theoretical argument which states 
that the high-Q 2 behavior of the three form factors 
should be the same. From the observations (i), (ii) (iii) 
one can already conclude that this property is indeed 
seen in the data if for F V1 the Belle data points are used. 
However, for a more quantitative analyses, the slight log- 
arithmic rise indicated by the data should be taken into 
account. 

To describe the logarithmic rise, the simplest use of the 
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Fig. 1: Form factors Q 2 Fp^iP = V^v') vs Q 2: experimental 
data from Cello and Cleo [J,|3| (black dots), BaBar [j| (red 
squares), and the data borrowed from the time- like region 
(green triangles). Dashed lines - the results from [3] which 
obey the factorization theorem at Q 2 — > oo; solid lines - our 
fits for 4 /=0) = r s = 0.05 GeV 2 . 



duality concept — a replacement of the absorptive part 
of the form factor beyond the resonance region by the 
lowest-order perturbative QCD spectral density — cannot 
be maintained: To account in a phenomenological way for 
higher-order effects, the lowest-order perturbative spec- 
tral density will be be multiplied by a correction factor 
R(s) which goes to one for high values of s. Here s de- 
notes the square of the energy variable. It turns out that 
it is sufficient to describe R(s) as a simple function of a 
single fit parameter r: R(s) = 1 — -. The function R(s) 
starts at the effective threshold relevant for each process. 
(For technical details see the next Section). 1 

Our simple model for the full set of the form factors 
thus involves six parameters: the three effective thresh- 
olds Sq 1 ^ , s s for the 7 = 0,7 = 1 and ss components 

of the form factors and the three r parameters rq~° , 
r s . The decay constants f s and f q and the r\ — 77' mixing 
angle 4> are taken from Ref. [26J . 



1 In principle, r could be Q 2 -dependent, e.g., r(Q 2 ) = 
This choice would restore saturation at large Q 2 . 
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Fig. 2: Form factor Q 2 F 7TJ (Q 2 ) vs Q 2 : experimental data 
from Cleo [2| (black dots), BaBar |^| (red squares) and Belle 
pl| (blue triangles); dashed line - the results from [l4| which 
obey the factorization theorem at Q 2 — > 00. Solid lines - our 
fits. In the fit to the Belle data, the value of r^ 7-1 ' is taken 
identical to r^ -0 ' for the nonstrange component of i ? r?7 and 
F v ' 1 in Fig. [T] The grey shaded region corresponds to the 

range 0.05 GeV 2 < r^ 1 ' < 0.17 GeV 2 . 



The form factors of 77 and r\' 

Let us first set r s = rq 1-0 ^ — 0, as implied by the 
pQCD factorization theorem. The fit to the existing data 
yields s q I=1) = 0.67±0.07 GeV 2 and s s = 1.0±0.03 GeV 2 
with x 2 /DOF = 72/71. Strictly speaking, in view of 
the experimental data, further improvements are unnec- 
essary. Nevertheless, having in mind the BaBar result for 
the pion form factor, let us allow also nonzero values of 
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r setting r q l °' = r s (at the present accuracy of the data 
it makes no sense to treat them independently). Then 

the fitting procedure gives Sq /_0 ^ = 0.57 ± 0.07 GeV 2 , 



s s = 0.94 ± 0.03 GeV 2 and r s = r q ' = 0.05 ± 0.01 
GeV 2 with \ 2 /DOF = 60/70. The corresponding results 
are shown in Fig. Q] Notice that the threshold values ob- 
tained by the fits are not far from Sq 1-1 ^ — 0.56 GeV 2 
and s s = 0.76 GeV 2 suggested from the relevant elements 
of the rj — rj' mass matrix [26j . 



The 7T° form factor 

A fit to the data set containing the Cleo and the 
BaBar results in Fig. pleads to a much larger paramater 
r ( /=1 ) = 0.14±0.014 GeV 2 than obtained for r q l=0) found 
above. The x value is \ 2 /DOF = 26/30 and the thresh- 
old parameter s q I=1) = 0.32 ± 0.07 GeV 2 . Most impor- 
tant, by setting A I=1 ^> = one gets an extremely bad 
fit with X 2 /DOF = 88/31. Obviously, the BaBar high- 
er 2 data are not compatible with the asymptotic pQCD 
result at large Q 2 (15-35 GeV 2 ). 

In contrast, a fit to the data set containing the Cleo and 
the Belle data in Fig. [2] is fully compatible with r, 



(i=i) 



0: for s 



(7=1) 



0.64 ±0.01 GeV 2 , one gets an excellent fit 



to the data with \ 2 = 16/29. No further improvements 
are necessary. Nevertheless, allowing in addition for a 



nonzero parameter r, 
GeV 2 and s q I=1) 



(i=i) 



leads to r. 



(i=i) 



0.06 ± 0.02 



0.5 ± 0.06 GeV 2 with x 2 /DOF = 
10/28. Remarkably, this value of rq 1-1 ^ is equal to rq /_0 ^ 
as obtained from the rj and rj' data. Also the effective 
thresholds for the nonstrange quark sector are very close 
to each other: Sq /_1 ' ~ s^ -0 " 1 . 

The shaded region in Fig. [2] corresponds to the varia- 
tion of r in the range 0.05 GeV 2 < r < 0.17 GeV 2 . 

To summarize, using the Belle data for the large-Q 2 
region, the previously puzzling difference between the 
(i],r)')j and ir-f form factors is no more present: all three 
processes — after taking particle mixing into account — 
can be well described by only two effective thresholds 
and a small universal parameter r. 

Certainly, more precise measurements are needed to 
establish a small logarithmic increase of Q 2 Fp 1 (Q 2 ) for 
high Q 2 values indicated by the data and parametrized 
by a nonzero r ~ 0.05 GeV 2 . 



general decomposition contains four independent Lorentz 
structures, but for our purpose only one structure is 
needed fl4j ] 



(2) 



The invariant amplitude F(p 2 ,Q 2 ) satisfies the spectral 
representations in p 2 at fixed Q 2 : 



oo 

F( P 2 ,Q 2 ) = - [ -^A( S ,Q 2 ), 

7T J S — £T 



(3) 



where A(s, Q 2 ) is the physical spectral density and sth 
denotes the physical threshold. Perturbation theory 
yields the spectral density as an expansion in powers of 
a s . The lowest order contribution, ApQ CD (s,Q 2 ), corre- 
sponds to the one-loop triangle diagram with the axial 
current and two vector currents in the vertices [20 22] 
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yl — 4m 2 /s. 



(4) 



Here m denotes the mass of the quark propagating in 
the loop. The integral of ApQ CD (s, Q 2 ) from s = 4m 2 to 
infinity is independent of m and Q 2 and gives the axial 
anomaly [23| 



ds A 



(°) ( s O 2 ) - — 



(5) 
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According to the Adler-Bardeen theorem [241 ] . radiative 
corrections to the anomaly vanish: Higher order QCD 
calculation can change A p Q CD but not the integral. 

Non-perturbative QCD effects strongly distort 
A(s, Q 2 ) compared with ApQ CD (s,<3 2 ) in the low-s 
region: A meson pole and the hadron continuum are 
generated. Nevertheless, the integral of the entire 
absorptive part A(s, Q 2 ) remains unchanged, still 
representing the anomaly: 



I dsA( s ,Q 2 ) = ±-. 



(6) 



TECHNICAL DETAILS 



We now provide some details of our calculation of the 
P7 form factors. For subtleties, we refer to (l7| . 
Our starting point is the amplitude 



(0|^ |7(<72)7* (<7i)> = e 2 T mP {p\ qi ,q 2 )e1el 

p = qi + <?2- 



(i) 



Here £1.2 denote the photon polarization vectors. This 
amplitude is considered for q 2 — —Q 2 and q\ — 0. Its 



For the case of the isovector uu ^ d axial current the spec- 
trum contains the 7r°-meson pole. The physical absorp- 
tive part of F(p 2 , Q 2 ) reads 

A(s,Q 2 )= _ (7) 
Trd(s - ml) V2U F^iQ 2 ) + 6(s - s th ) A^(s, Q 2 ). 

Here A^^s, Q 2 ) denotes the hadron-continuum contri- 
bution in the isovector channel. In ((7|), the 7r7 form factor 
we are interested in appears together with the 7r meson 
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(5-function and the pion decay constant. The anomaly 
sum rule for F 7T ^(Q 2 ) then takes the form 



F^(Q 2 ) = 



2V2 7T 2 U 



1 - 2vr 



ds A^(a,Q 2 ) 



(8) 



For the 777 and ^'7 form factors, one has to consider the 
isoscalar currents qq = (uu + dd)/y/2 and ss, separately. 
The formulae for F qq {Q 2 ) and F Ss (Q 2 ) are identical to 
(jHJ except for the replacements v^/^ by f q and / s and 
A c™t b y A lont and A fo„t> respectively. 

For each channel, the relevant threshold Sth should be 
used. Taking 77 — 77' mixing [25|, [26| into account leads to 

F V1 (Q 2 ) = -^=F qq (Q 2 ) cos0-^ Ss (Q 2 ) sin0, 

P/ 7 (Q 2 ) = ~^=F qq (Q 2 ) smq>+ ±F SS (Q 2 ) cos^. (9) 

The 77 — 77' mixing angle is known to be (f> ~ 39°; the 
decay constants are taken to be f q = 1.07 f n , f s — 1.36/^ 

M- 

According to ([8]) and ((9]), the calculation of the P7 
form factors requires an Ansatz for the continuum spec- 
tral densities A cont (s,<3 2 ) for all three cases. 

The quark-hadron duality suggests that at large val- 
ues of s, above the resonance region, the hadron spectral 
density should be very close to the perturbative QCD 
spectral density. We therefore use the simple Ansatz 



A con t(s, Q 2 
R(s o 



(s-s th )R(s)A§l u (s,Q 



1. 



(10) 



It turns out [17j that for the large-Q 2 behavior of the 
form factor the behaviour of R(s) at large s is essential: 
e.g. in order to have the logarithmic rise of Q 2 F(Q 2 ), 
R(s) should contain a 1/s-correction: For R(s) = 1 — r/s 
starting at a finite energy sq, Sq > s t h, one finds 



Q 2 F(Q 2 ) 



Q 2 



Q 2 
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r log 
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(11) 



Notice that the part of the integral © from s t h to so 
scales as 1/Q 2 and is therefore not relevant for the effect 
discussed here. Thus, neither the details of R(s) at small 
s nor the presence of higher powers of 1/s affect this 



behavior. The negative linear term — r/s in R(s) is re- 
sponsible for the logarithmic increase. If absent, the form 
factor scales as F(Q 2 ) ~ 1/Q 2 - Evidently, for a universal 
r the high-Q 2 behavior of all three form factors turns out 
to be the same. For small r, the differences due to differ- 
ent thresholds and decay constants affect the magnitude 
of the form factors but not their slope. 

3. CONCLUSIONS 

We revisited the P7 transition form factors, P = 
77,77,77', with special emphasis on the new data on the 
7T7 form factor reported by Belle Use is made of 

the exact anomaly sum rule which relates the integral 
over the hadron spectrum to the axial anomaly. This 
approach has the advantage that the QCD factorization 
theorem and the meson distribution amplitudes do not 
enter the analyses. Thus, the three processes can be eas- 
ily compared with each other. 

• We report that the 77* — > P form factors of 7T, 77 and 77' 
are fully compatible with each other — if for the ttj form 
factor the recent Belle data are applied: the parameter r 
and the effective threshold used in the description for the 
1 = nonstrange continuum and the ones for the 1 = 1 
continuum agree with each other. Thus, the Belle data 
resolve the puzzle of a qualitatively different behavior of 
the nonstrange component in 77 and rf on one hand and 
of 7T on the other hand. 

• The Belle data for the 777 form factor are compatible 
with the asymptotic pQCD formula indicating that cor- 
rections to the asymptotic behavior are small already at 
Q 2 > 10 GeV 2 . 

• Still, our best fits to the data — for all three processes — 
suggest a slight increase of the product Q 2 F{Q 2 ) at high 
Q 2 . If confirmed by future experiments, this would put 
QCD factorization into question and would suggest that 
the full spectral density of the dispersion representation 
for the form factor is dual to the lowest order pQCD spec- 
tral density only by including an effective 1 / s-correction 
term. 

Acknowledgments. We are grateful to W. Lucha, 
B. Moussallam, J. Pawlowski, H. Sazdjian, and 
O. Teryaev for valuable discussions. D. M. was supported 
by the Austrian Science Fund (FWF) under Project 
No. P22843. 



[1] CELLO Collaboration, H. J. Behrend et al, 

Z. Phys. C49, 401 (1991). 
[2] CLEO Collaboration, J. Gronberg et al., Phys. Rev. D57, 

33 (1998). 

[3] BaBar Collaboration, B. Aubert et al., Phys. Rev. D74, 
012002 (2006). 

[4] BaBar Collaboration, B. Aubert et al., Phys. Rev. D80, 



052002 (2009). 

[5] BaBar Collaboration, P. del Amo Sanchez, 

Phys. Rev. D84, 052001 (2011). 
[6] A. V. Radyushkin, Phys. Rev. D80, 094009 (2009). 
[7] H. L. L. Roberts, C. D. Roberts, A. Bashir, 

L. X. Gutierrez-Guerrero, P. C. Tandy, Phys. Rev. C82, 

065202 (2010). 



5 



[8] A. Dorokhov, JETP Lett. 91, 163 (2010). 

[9] S. S. Agaev, V. M. Braun, N. Offen , and F. A. Po rkert, 

Phys. Rev. D83, 054020 (2011); larXiv: 1206.39681 [hep- 

ph]. 

[10] Y. N. Klopot, A. G. Oganesian, and O. V. Teryaev, 
Phys. Lett. B695, 130 (2011); Phys. Rev. D84, 051901 
(2011). 

[11] S. J. Brodsky, F.-G. Cao, and G. F. de Teramond, 
Phys. Rev. D84, 033001 (2011); Phys. Rev. D84, 075012 
(2011). 

[12] P. Kroll, Eur. Phys. J. C71 1623 (2011). 

[13] A. P. Bakulev, S. V. Mikhailov, A. V. Pimikov, 

and N. G. Stefanis, Phys. Rev. D84, 034014 (2011); 

Phys. Rev. D86, 031501 (2012). 
[14] I. Balakireva, W. Lucha, and D. Melikhov, 

Phys. Rev. D85, 036006 (2 012); J. Phys. G39, 

055007 (2012) |arXiv:1103.378 1 ; W. Lu cha and D. Me - 

likhov, J. Phys. G39, 045003 (201 2) [arX"iv:1110.2080l ; 

Phys. Rev. D86, 016001 (2012) arXiv: 1205.4587^. 
[15] H. Czyz, S. Ivashyn, A. Korchin, and O. Shekhovtsova, 

Phys. Rev. D85 094010 (2012). 
[16] C.-C. Lih and C.-Q. Ceng, Phys. Rev. C85, 018201 



(2012). 

[17] D. Melikhov and B. Stech, Phys. Rev. D85, 051901 
(2012). 

[18] Belle Collaboration, S. Uehara et al, arXiv:1205.3249l 
[hep- ex] 

[19] G. P. Lepage and S. J. Brodsky, Phys. Rev. D22, 2157 
(1980). 

[20] J. Horejsi and O. Teryaev, Z. Phys. C 65, 691 (1995). 
[21] D. Melikhov and B. Stech, Phys. Rev. Lett. 88, 151601 
(2002). 

[22] D. Melikhov, Phys. Lett. B 380 , 363 (1996); 

Eur. Phys. J. direct C 4, 2 (2002) [hep-ph/0110087] . 
[23] S. Adler, Phys. Rev. 177, 2426 (1969); J. S. Bell and 

R. Jackiw, Nuovo Cimento 60A, 47 (1969). 
[24] S. Adler and B. Bardeen, Phys. Rev. 182, 1517 (1969). 
[25] V. V. Anisovich, D. I. Melikhov, and V. A. Nikonov, 

Phys. Rev. D55, 2918 (1997); V. V. Anisovich, 

D. V. Bugg, D. I. Melikhov, V. A. Nikonov, Phys. Lett. B 

404, 166 (1997). 
[26] T. Feldmann, P. Kroll, and B. Stech, Phys. Rev. D58, 

114006 (1998); Phys. Lett. B449, 339 (1999). 



